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Summary 


The research during the second reporting period has focused on robust control design for 
hypersonic vehicles. An already existing design for the Hypersonic Winged-Cone Con- 
figuration has been enhanced. Uncertainty models for the effects of propulsion system 
perturbations due to angle of attack variations, structural vibrations, and uncertainty in 
control effectiveness were developed. Using Hoo and //-synthesis techniques, various con- 
trol designs were performed in order to investigate the impact of these effects on achievable 
robust performance. 
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1 Introduction 


The objective of this research is to address the issues associated with the design of ro- 
bust integrated flight control systems for future hypersonic vehicles with airbreathing 
propulsion systems. It is anticipated that such vehicles will exhibit significant interac- 
tions between rigid body (airframe) dynamics, structural dynamics and engine dynamics. 
The uncertainty in the initial dynamic models developed for these vehicles will also be 
high. The main reason that highly interactive uncertain dynamics are to be expected is 
that scramjet engines will be the primary source of propulsion at hypervelocity speeds. 
Wind-tunnel testing as a result will be limited, and it will be necessary to gain experience 
in actual flight testing of such vehicles. This means that initial flight control system design 
efforts will rely more heavily on theoretical and computer based models, than has been 
the case for subsonic and supersonic aircraft. Also, propulsion system sensitivity to angle 
of attack variations and structural vibrations will lead to highly interactive dynamics. 


In this study, the current major research issues from a flight control viewpoint are: (1) the 
development of models that are representative of the interactive dynamics that can occur 
in such vehicles, (2) the development of representative uncertainty models for these dy- 
namics and (3) the development of practical approaches to designing multivariable flight 
control systems that guarantee adequate performance in the presence of uncertainty. The 
research done during the second reporting period has been focusing on items (2) and (3). 


The hypersonic vehicle model used in this study [1] neglects both the effects of angle of 
attack variations on propulsion system performance and of elastic body bending. For long, 
slender bodies with a considerable amount of the compression of the flow going through 
the engine taking place on the forebody, this assumption may not necessarily be valid 
for the entire flight regime. Changing the angle of attack for control purposes alters the 
forebody flowfield. This effect propagates through the engine and results in variations in 
thrust vector magnitude and direction. Low structural vibrations frequencies may lead 
to considerable elastic-rigid body mode interactions while flexible body bending of fore- 
and aftbody again affects the flowfield relevant for propulsion system performance. 


In the framework of robust control design these effects are treated as uncertainties. This 
study has resulted in uncertainty models capturing the individual characteristics of these 
phenomena. The models were developed and incorporated into a control design structure 
which evolved from an earlier study described in Ref. [2]. A variety of robust controllers 
were designed utilizing Hoo and ^/-synthesis techniques and the sensitivity of achievable 
robust performance to the introduced uncertainty levels were investigated. A thorough 
description as well as a comprehensive discussion of the results was presented in Ref. [3] 
at this year’s AIAA Guidance, Navigation and Control Conference in Monterey, CA, (see 
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Appendix). Therefore, the emphasis in this report is on a comparison of our design to 
the design in Ref. [2], and on additional background information not given in the paper. 


2 Framework for the Nominal Design 

The vehicle model used in this study is the Winged-Cone Configuration described in Ref. 
[1]. To perform control studies a five state linear model was provided by NASA Langley 
Research Center (courtesy of J.D. McMinn). This model was derived at an accelerated 
flight condition at Mach 8 with velocity, angle of attack, pitch rate, pitch attitude and 
altitude as state variables, and with elevon deflection and fuel equivalence ratio as control 
variables. The outline of our controller design structure is based on a problem formulation 
given in Ref. [2]. The modifications and enhancements are described in the following. 


Under the assumption that normal acceleration measurements using accelerometers are 
easier to obtain than angle of attack measurements, a is replaced by the normal acceler- 
ation n z as a measurement signal which is fed back to the controller. 


The time responses in angle of attack and control deflections given in Ref. [2] for an Hoc- 
design encountering longitudinal and vertical turbulence exhibit a rather noisy behavior. 
In order to obtain an improved turbulence attenuation, a variety of Hoo-designs with dif- 
ferent weights (constant as well as frequency dependent) on angle of attack, pitch rate, 
and normal acceleration were examined. This approach turned out to be ineffective in 
attenuating the influence of atmospheric turbulence. It was determined that an increased 
penalty on the control rates 8e and 8t) is effective in limiting the turbulence impact on 
the flight behavior of the vehicle. Accordingly, the corresponding weights were increased 
to 


w it = 90 (1) 

l Vi, = 50 (2) 

with respect to the weights as defined in Ref. [2]. Also in our study (see Ap- 
pendix) all weights on the performance outputs are scaled in order to obtain an oo- 
norm from disturbance inputs to performance outputs for a suboptimal Hoo-controller 
of || T e d ||oo= 0.78 < 1. Moreover, angle of attack, pitch rate and pitch attitude were 
removed as performance variables. 
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The effect of turbulence on the system is modeled by 


where 


x = Ax + Bu + B g 

1 

1 




9s J 


—a n 

— 012 

0 


— O21 

— <*22 

0 

b 9 = 

®31 

0 

— a 32 
0 

—<*33 

0 


0 

0 

0 


( 3 ) 


( 4 ) 


The elements in B g are the corresponding elements of the A matrix. The gust components 
in angle of attack and pitch rate are 
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( 5 ) 

( 6 ) 


respectively. This turbulence model is taken from Ref. [4]. The gust terms V g , w g and 
w g are provided by a Dryden turbulence model. This differs from the modeling used in 
Ref. [2] in that the elements a 32 and a 33 in Eq. 4 are zero in the formulation in Ref. [2]. 
These terms represent the influence of a g and q g in the pitching moment. 


To illustrate the differences in the two designs a comparison of the singular values of the 
loop transfer function with the loop broken where velocity and altitude are fed back to 
the controller is shown in Fig. 1. In both cases Hoo-controllers for the nominal model 
without uncertainty were used. It can be seen that the bandwidth of our model is reduced 
compared to the design in Ref. [2]. 


3 Uncertainty Modeling 

The main task during the second reporting period was the development of uncertainty 
models representing propulsive and aeroelastic effects. A detailed description of what 
these models are and how they are implemented in the controller design is given in the 
Appendix. Some of the motivational aspects and background information is addressed 
below. 
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Figure 1: Comparison of singular value plots for loop broken at the V and h outputs, 
design in Ref. [2] (dash) vs. our design (solid). 

Angle of attack changes due to vehicle control alters the flowfield along the forebody of 
a hypersonic vehicle. These variations in the flowfield properties propagate through the 
scramjet engines and affect the propulsion system performance. This results in perturba- 
tions on the thrust vector acting on the vehicle and thus influences stability and control 
[5]. During our first reporting period the HYTHRUST code was used to address this 
issue. It was observed that for the Winged-Cone Configuration the most sensitive vari- 
able to these propulsive variations was the pitching moment. At certain Mach numbers 
the contribution of the propulsion system to the change in the overall pitching moment 
was of the same order of magnitude as the aerodynamic contribution of the basic vehicle 
for a given angle of attack change. Hence, the pitching moment sensitivity to angle of 
attack variations cm a = dcM/da is considered the variable which will be affected the 
most by propulsion system perturbations. In robust control design, this effect is treated 
as parameter uncertainty in CMa- 


Aeroelastic effects will also be of great significance when designing a flight control system 
for hypersonic vehicles. Elevon deflections will excite flexible modes of the configura- 
tion resulting in elastic deformations of the vehicle. Bending of fore- and aftbody again 
changes the flowfield properties and inflow conditions of the propulsion system leading to 
the effects described above. Since no aeroelastic information of the Winged-Cone Vehicle 
was available, a simple yet complete technique was developed to model the uncertainty 
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due to this effect. 


The long, slender shape of the hypersonic configuration was approximated by a 2- 
dimensional uniform beam. The flexural deflections of such a beam are governed by 
the partial differential equation 


El = p{x, t) 


dx 4 


dt 2 


( 7 ) 


The longitudinal coordinate is x , the vertical coordinate (deflection) is y — /(x,t), t is 
time, E is the Young’s modulus of elasticity, I is the moment of inertia, m is the mass per 
unit length, and p is the load per unit length and time. Using the principle of separation 
of variables, the general solution can be expressed by the summation of the products of 
the normal modes <j>i(x) multiplied by the factors Zi(t ) 


y(x,t) = Y^<f>i(x)zi(t) i= 1,2,.... (8) 

i 

The normal modes <f>i(x) are characterized by free vibration (p(x,t) = 0) and satisfy the 
differential equation 


El <f>iM(x) Zi{t) _ 

m <^,(x) Zi(t) ' 

( 9 ) 

which is obtained by applying the principle of separation of variables on Eq. 7. Substi- 
tuting Eqs. 8 and 9 into Eq. 7, integrating over the length of the beam L and observing 
the orthogonality property of normal modes ( <f> n • <f> m = 0 if m) yields 

u>?Zi(t) f m<f) 2 i {x)dx— f <f>i(x) p(x,t) dx — 5,-(f) f m<j) l 2 (x)dx . 
Jo Jo Jo 

(10) 

Defining modal mass as 


M i = / m <f> 2 (x) dx 
Jo 

(11) 

and modal force as 


Fi = f <f>i p[x yt'j dx , 

(12) 

Eq. 10 can be written as 



Mi Zi{t) + us] M, zi{t) = Fi{t ) , i = 1,2, . . . (13) 

Introducing modal damping by simply adding the damping term gives 

F (A 

z(t) + 2 (iUiZi(t) + uf Zi(t) = . (14) 
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Since the elastic vibrations are caused by elevon deflections the force acting on the body 
is the lift increment at the elevon location: 


p(x,t ) = P(t ) at x = x p . 


So the modal force simplifies to 


Fi(t) = P(t ) <f>i{x p ) . 


(15) 

(16) 


If we further replace the deflection y(x,t) by the deflection angle 6(x y t) using 


6(x,t) ^ tan 6{x,t) = (17) 

Ox 

and the principle of seperation of variables, we obtain 

*(*,<) = £[*(*> *iw] = • < 18 > 

Substituting Eqs. 16 and 18 into Eq. 14 and applying the Laplace transform yields a 
second order transfer function from the elevon lift force P(s) to the body deflection angle 
8(x, s ) in terms of the elastic characteristics of the body: 


Oi(x,s) 


H x p) 


d<t>,(x) 

dx 


(s 2 + 2 £,-il>,s + Mi 


P(s) 


(19) 


where i accounts for the individual modes. The implementation of this transfer function 
into the controller design as an uncertainty model is described in the Appendix. 


The natural frequencies of a hypersonic configuration will most likely be uncertain to 
begin with. Additionally, aerodynamic heating during high speed atmospheric flight will 
affect the elastic characteristics of the vehicle. Increased heating loads tend to decrease 
the natural frequencies of the flexible modes. The transfer function in Eq. 19 possesses 
the property that for decreasing natural frequencies, the peak value of its frequency re- 
sponse increases. The frequency response for two different sets of natural frequencies is 
shown in Fig. 2. As the elastic mode frequencies approach the rigid body frequencies, the 
level of uncertainty increases. 


Uncertainty in the control effectiveness was already addressed in Ref. [2]. There, con- 
stant levels of uncertainty were inserted simultaneously in both control channels. In our 
study we extend this approach to frequency dependent uncertainty models to account 
for increasing uncertainty with increasing frequency at the plant input. The uncertainty 
weight was chosen to be 

TX7 &act(s+10) t ^ 

Wi “' = , + 1000 (20) 


where k act determines the uncertainty level for low frequencies. This weighting function 
is illustrated in Fig. 3. 
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4 Results and Conclusions 


A complete discussion of the robust control design procedure and the results is given in 
the Appendix. The main conclusions are summarized below: 


• The design for robust performance with all uncertainties modeled simultaneously 
demonstrates that Hoc and // analysis and synthesis methods are ideal for addressing 
multiple uncertainty sources in hypersonic flight control design. 

• Flexible fuselage bending and control effectiveness has a greater impact on achiev- 
able robust performance in comparison to uncertainty due to propulsion system 
variations. To a certain extend, this confirms the observations made with the 
HYTHRUST code during our first reporting period [5]. 

• Mutual coupling effects among the individual phenomena decrease the admissible 
uncertainty levels for which robust performance is achievable. This shows that 
sophisticated uncertainty modeling with reduced conservatism will be crucial for a 
successful flight control design for hypersonic vehicles. 

• The controllers obtained using //-synthesis techniques are usually of rather high or- 
der. An attempt to reduce the order of the controller using various model reduction 
techniques failed in that robust performance was lost for the system using the re- 
duced order controllers. This shows the need to employ fixed order controller design 
techniques. 

5 Future Research 

Based on the results from our first year’s research the focus of our future work will be on 
three major topics: 

1. Fixed order controller design methodology will be developed and applied on the 
hypersonic vehicle model investigated in this report. First steps have already been 
taken to extend the fixed order Hoo design described in Ref. [6] to a //-synthesis 
framework and encouraging results were obtained for a simple example. 

2. Atmospheric turbulence has a significant impact on the level of control activity in 
hypersonic flight. A possibility to attenuate this influence is to impose an additional 
H 2 constraint. This will be done in the differential game setting for fixed order Hoo 
control described in Refs. [6] and [7]. This is the so-called mixed H 2 /Hoo problem 
for which no complete solution presently exists even in the full order case. 
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3. To further reduce conservatism in the uncertainty modeling, the issue of real pa- 
rameter uncertainty will be investigated, //-synthesis accounts for structure in the 
uncertainty but treats real parameters as complex which is overly conservative. 
Again, the differential game formulation can be employed to treat uncertain pa- 
rameter variations (in addition to external disturbances) as an opponent trying to 
maximize the performance index while the controller tries to minimize the same 
performance index. Another approach to solve the real // problem that is currently 
under development is the mixed //-analysis using the so-called G-scales in addition 
to the D-scales which is described in Refs. [8] and [9]. Also, results have already 
been obtained by treating the uncertainty as nonlinearity and applying techniques 
from Absolute Stability Theory such as the Popov Criterion, e.g. see Ref. [10]. We 
will look into these different approaches as well as possible combinations in order to 
solve the problem of real/mixed //-synthesis. 
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Abstract 

The influence of propulsion system variations and 
elastic fuselage behavior on the flight control system 
of an airbreathing hypersonic vehicle is investigated. 
Thrust vector magnitude and direction changes due 
to angle of attack variations affect the pitching mo- 
ment. Low structural vibration frequencies may occur 
close to the rigid body modes influencing the angle of 
attack and lead to possible cross coupling. These ef- 
fects are modeled as uncertainties in the context of 
a robust control study of a hypersonic vehicle model 
accelerating through Mach 8 using Hoc cm d \x syn- 
thesis techniques. Various levels of uncertainty are 
introduced into the system. Both individual and si- 
multaneous appearance of uncertainty are considered. 
The results indicate that the chosen design technique 
is suitable for this kind of problem provided that a 
fairly good knowledge of the effects mentioned above 
is available. The order of the designed controller is 
reduced but robust performance is lost which shows 
the need for fixed order design techniques. 

1. Introduction 

Hypersonic atmospheric flight will be one of the 
most challenging efforts undertaken by aerospace sci- 
entists at the end of this century. Aerospace vehicles 
have to perform in speed ranges from subsonic to hy- 
personic up to Mach 25, and therefore will encounter 
a variety of effects that have not been addressed be- 
fore in the development of one single vehicle. One of 
the most crucial aspects will be a strong interaction 
between aerodynamics, structure, and propulsion sys- 
tem and their impact on performance, guidance and 
control characteristics [1]. 

In hypersonic propulsion, unlike present sub- 
sonic/supersonic engines, the entire vehicle aerody- 
namic configuration must be considered as part of the 
propulsion system. Before reaching the inlet the ex- 
ternal flow will be compressed along the forebody, and 
after leaving the combustor the aftbody is utilized as 
an external nozzle. Changing the angle of attack for 
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control purposes alters these flow fields and results in 
variations in thrust vector magnitude and direction 
influencing stability and control of the vehicle. The 
most sensitive quantity in this respect is the pitching 
moment [2], the control of which is vital to stabilize 
typically long and slender bodies of aerospace vehicles 
at this speed regime. 

The need for a low structural weight mass fraction 
introduces relatively low structural vibration frequen- 
cies. Significant elastic-rigid body mode interactions 
are likely to occur imposing additional requirements 
on the flight control system. Additionally, bending of 
the forebody influences the flow conditions at the in- 
let which propagates through the engine and together 
with the elastic deformation of the aftbody again af- 
fects the thrust vector [3]. Determination of the elas- 
tic mode shapes will be highly uncertain. 

These features, among others, illustrate the neces- 
sity of designing a robust, highly integrated flight con- 
trol system for future hypersonic vehicles with air- 
breathing propulsion systems. In this paper, uncer- 
tainty models are developed to address the two ma- 
jor effects mentioned above, and, together with addi- 
tional uncertainty on the control effectiveness, their 
impact on achievable robust perform since is investi- 
gated. This is done within the context of a flight 
control system design using /j-synthesis design tech- 
niques. The uncertainty modeling is described, and 
a sensitivity study of achievable robust performance 
to the uncertainty level is carried out. In an addi- 
tional design step the order of the obtained controller 
is reduced using several different methods, and a com- 
parison with the full order design is presented. 

2. Brief Review of /i-Syn thesis 

This section provides a concise summary of the gen- 
eral framework of /i-synthesis in order to define the 
terminology used in the paper. A thorough discus- 
sion on Hoo and /i control theory is provided in Refs. 
[4]-[6]- 

The framework for robust control design is illus- 
trated in Fig. 1. P is the generalized plant, K 
denotes the controller, and A represents the uncer- 
tainties on the system due to parameter variations, 
unmodeled dynamics, nonlinearities and various as- 
sumptions. The input vector, d, denotes all external 
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Figure 1: Framework for robust control design. 

disturbances acting on the plant, such as commands, 
disturbances and sensor noise. The output vector e 
contains the regulated performance variables. Feed- 
back measurements are denoted by the vector y, and 
« represents the control vector that is the input to 
the plant actuators. The inputs and outputs corre- 
sponding to the uncertainties in the system are w and 
z y respectively. The input disturbances and uncer- 
tainties are normalized to unity bounds by constant 
or frequency dependent weights, which are incorpo- 
rated into the generalized plant, P. The uncertainty 
A is considered to be complex with a block diagonal 
structure consisting of s scalar blocks and / full blocks 

A = diagfMn • ..6,I rt ,Ai . .A/} (1) 

whf>rp 

6, € C, Aj € C mjxm> . (2) 

The bounded set A < is defined as 

A, = {A: |6<| < 6, &(Aj) < 5} 

= {A:*(A)<6} (3) 

If the controller is wrapped into the plant P, the re- 
sulting system M is defined by the lower Linear Frac- 
tional Transformation 

M = Fi(PyK) 

= Pn + PnK(I-P 2 2K)- 1 Pn (4) 

This defines the analysis structure which is shown in 
Fig. 2. It is assumed that M is stable. 

To account for the structure in the uncertainty A, 
the Structured Singular Value fi(Mn) is introduced 
which is defined by 

{ {min A€ Aoo [ff(A): det (/ - M n A) = 0]} 

0 if det(/-MnA)^0,VA€ Aoo. 

( 5 ) 



Figure 2: Analysis structure. 

The corresponding //-measure is 

H-WiilU = sup/i(Mn) (6) 

U> 

Note .that the //-measure is not a norm and depends 
both on M and the structure of A. 

With these preliminaries stated the concepts of per- 
formance, stability, and robustness can now be ad- 
dressed. Referring to Fig. 2, in //-synthesis the per- 
formance output e is scaled so that nominal perfor- 
mance (NP) is said to be satisfied if the disturbance 
attenuation condition 

IIAMU < 1 (7) 

is satisfied. In Eq. (7) the Hoo-norm is defined as 

IIM 22 II 00 = sup a (M 2 2 (»w)) . (8) 

W 

Robust stability (RS) is satisfied if and only if (iff) 

||M u |U<l,Ae Al (9) 

For unstructured uncertainty, the //-measure equals 
the Hoc -norm. 

The robust performance problem is converted into 
an equivalent robust stability problem by augmenting 
the uncertainty block with a fictitious full uncertainty 
block A p as shown in Fig. 3. With this formulation, 
robust performance (RP) is satisfied iff 

II^IU < 1 * Aouj € Aj . (10) 

where A a u s stands for the augmented uncertainty. 

The Structured Singular Value cannot be calculated 
directly. An upper bound can be* computed from: 

/x(Af) < inf a(DMD~ l ) (11) 

i 

where D = diag{d,/j} with the same dimensional 
structure as A. The equality on this upper bound 
holds in some cases, but is not true in general. How- 
ever, the optimization problem for the upper bound is 
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Figure 3: Structure for addressing robust perfor- 
mance. 

convex leading to an upper bound value that is fairly 
close to the /i value in most cases. The robust per- 
formance test can be expressed in terms of this upper 
bound: 

inf \\DMD~ 1 ||oo < 1 =* RP. (12) 

Note that it has become an oo-norm test. 

Hoo-desigB and /i- analysis are combined to perform 
/i-synthesis. Recalling Eq. (4), an optimal controller 
can be achieved by solving the optimization problem 

mf ||DF,(P,tf)D- J ||=o (13) 

Jv, D 

Presently, it is not known how to solve this minimiza- 
tion. An approximation to /i-synthesis involves a se- 
quence of minimizations, first over the controller vari- 
able K (holding the D variable fixed), and then over 
the D variable (holding the K variable fixed). This is 
often referred to as the D-K iteration. 

3. Hypersonic Vehicle Model 


Since the linear model is obtained at a non- 
equilibrium flight condition where the vehicle accel- 
eration is non-zero, it should be noted that the state 
and control variables are perturbation quantities rep- 
resenting deviations from a climbing and accelerat- 
ing flight condition. This model represents pure rigid 
body dynamics and therefore does not account for any 
aeroelastic effects. Also, the propulsion system model 
used for the Winged-Cone Configuration does not in- 
clude sensitivity to angle of attack variations. On the 
other hand, these effects are considered crucial for ve- 
hicle stability and control and will be treated in this 
paper as uncertainty. 

The interconnection structure for the controller de- 
sign used in this study is shown in Fig. 4 and has been 
derived from the design in Ref. [8]. Velocity, altitude, 
normal acceleration, pitch rate, and pitch attitude are 
the fedback variables. Normal acceleration is deter- 
mined by 


= 5 ^ (i ‘ q) (16) 

where Vo is the velocity at the flight condition for 
which the linear model was obtained, 57.3 is the con- 
version from radians to degrees, and g is acceleration 
due to gravity. Control actuator dynamics are repre- 
sented by first order filters with 30 rad/sec bandwidth 
for elevon and 100 rad/sec bandwidth for fuel equiv- 
alence ratio. The design is carried out as a velocity 
and altitude command tracking system. Additional 
disturbances include Dryden turbulence models and 
sensor noise. 

The turbulence spectrum is defined by the weights 


*■-(«) 


Wort l 

L u s + Vo/L u 


(17) 


The hypersonic vehicle model used in this study is 
the Winged-Cone Configuration described in Ref. [7]. 
Main characteristics of this vehicle are an axisymmet- 
ric conical forebody, a cylindrical engine nacelle sec- 
tion with engine modules all around the body, and 
a cone frustrum engine nozzle section. In order to 
carry out control studies, a five state linear model of 
the longitudinal dynamics is used representing flight 
conditions for an accelerated flight through Mach 8 
at approximately 86000 ft (see Appendix). State and 
control variables are: 


■ V ' 


velocity (ft/sec) 

a 


angle of attack (deg) 

5 

= 

pitch rate (deg/sec) 

e 


pitch attitude (deg) 

h m 


altitude (ft) 


(14) 


[ 6e 1 _ [ symmetric elevon (deg) 
6n \ ~ [ fuel equivalence ratio (-) 


(15) 


for longitudinal and 


F„(«) 


j mrt 
V U 


fi + Vo/(y/3M 

(s + Vo/M 2 


(18) 


for vertical wind gusts [9]. For the reference altitude 
used in this study [10] 


t t u = 10.8 ft/sec, a w = 6.88 ft/sec, (19) 

L u = 65574 ft, L w = 26229 ft. (20) 

Turbulence introduces the gust quantities V g , a g and 
a g into the system (see Fig. 4 and Appendix). A 
measurement weight 


Wn<rf,e = 10- 6 I 5 (21) 

is imposed on the fedback variables. 

The weighting functions on velocity and altitude er- 
ror reflect specific performance requirements and were 
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Figure 4: Interconnection structure for controller design. 


derived in Ref. [8] from a near fuel optimum ascent 
trajectory for this vehicle. The weighting on the ve- 
locity and altitude errors are 


Wv = 

2.5 -10- 2 (s + 4.33 lO" 2 ) 

(22) 

* + 4.33- 10- 5 

w h = 

2.5 10- 2 (« + 4.95 10" 2 ) 

(23) 

s + 4.95 ■ 10 -5 


and allow for 10% overshoot and 5% steady state er- 
ror in the time responses of each velocity and altitude. 
Control deflections and control rates are the only 
other weighted performance outputs. The control de- 
flection weights are chosen to be 

W u = 1.25 (24) 

W* n = 3. (25) 

In order to reduce the sensitivity of the control re- 
sponse to atmospheric disturbances the weights on the 
actuator rates are selected to be 

W ie = 4.5 (26) 


W f „ = 2.5. (27) 

The uncertainties appearing at various locations in 
the diagram in Fig. 4 will be addressed in the follow- 
ing section. 

4. Uncertainty Modeling 

Aerospace vehicles in hypersonic flight regimes will 
typically utilize scramjet propulsion systems which 
are highly integrated into the airframe. This results in 
an increased sensitivity to variations in angle of attack 
[11]. The most important impact of these propulsive 
perturbations is on the pitching moment leading to 
significant control surface deflections to stabilize the 
vehicle [2]. Thi6 phenomenon is addressed as para- 
metric uncertainty in the pitching moment sensitivity 
to angle of attack variations, cm a- The uncertainty 
is represented by a scalar perturbation to the nom- 
inal model. This perturbation can be rearranged to 
obtain input and output quantities and A in 
this case is a scalar since only one parameter is per- 
turbed. The matrices C&, and D& can be chosen 
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from the uncertain model 


x = Ax + Ha + Bu 
z& = Ca* + £>a« 
y = Cx + Dti 
= A 2 a 


(28) 


If a scaling has been applied so that the values of A 
range between -1 and 1, Eq. (28) represents the set 
of all possible models. For the case of uncertainty in 
the parameter cm q which occurs in the (3,2) element 
of the A matrix, the uncertainty matrices are 


Ba 

= [001 0 0]^ 

(29) 

Ca 

= [01 000] 

(30) 

Da 

= 0 

(31) 


It should be noted that this approach leads to a con- 
servative representation since it also permits complex 
perturbations in real aerodynamic coefficients. This 
uncertainty is labeled “uncertainty 1” in Fig. 4. 

Significant coupling between the elastic and rigid 
body modes is expected for this vehicle type which 
has to be considereed in the design of a flight control 
system. Moreover, fuselage bending affects propul- 
sion system performance which in turn influences the 
rigid body flight dynamics or even excites the elas- 
tic modes [3]. Since the model used in this study 
comprises only rigid body dynamics and no aeroelas- 
tic information on the Winged-Cone Accelerator was 
available a simple yet concise method had to be de- 
veloped to introduce aeroelastic effects as uncertainty 
into the rigid body behavior. 

A second order transfer function was derived ex- 
hibiting elastic deformations of a long, slender, uni- 
form body caused by a force P(s) suddenly applied at 
a certain location of the body, x p [12]: 


( . # t (x) 

*<(*> s) = (fi2 + id ns + «*) Mi P{s) 


(32) 


0j(x,s) is the angle by which the deformed body is 
deflected from its original shape due to the I th elastic 
mode, ^i(x) is the mode shape function of the X th 
elastic mode, and u>i are the corresponding damping 
ratio and natural frequency, respectively, and Mi is 
the modal mass defined by 


Mi = 



where m(x) is the mass distribution and L the length 
of the body. In the case of the Winged Cone Vehicle 
structural excitement is caused by lift increments due 
to elevon deflections: 


P(s) = £L($) 

= 2 ^±qSre,6e(s) (34) 


where q is the dynamic pressure for the given flight 
condition, S r ef the reference area, and dci^a/ dSe the 
sensitivity of the lift coefficient of one elevon to the 
elevon deflection. This was graphically determined 
from Ref. [7], The body deflection angle 0(x,s) can 
also be interpreted as change in the angle of attack 
Aa(x,s) due to elastic deformation. As mentioned 
above, angle of attack perturbations affect the rigid 
body dynamics via the propulsion system, thus the 
transfer function given in Eq. (32) is evaluated at the 
location of the inlet entrance. The resulting estima- 
tion of the angle of attack variation due to flexible 
body motion excited by elevon deflections is then fed 
as uncertainty into the rigid body model, using the 
second column of the nominal plant A matrix (see 
Fig. 4 and Eq. (14)). 

Unfortunately, no information on the mode shapes 
of the Winged-Cone Accelerator was available 
for this study. So the rather crude assumption of 
a uniform beam with both ends free was chosen to 
model the elastic mode shapes. The analytical solu- 
tion to this problem can be found in any textbook on 
structural dynamics (e.g. [12]). The natural frequen- 
cies are taken from a generic NASP configuration in- 
vestigated in Ref. [13]. The first three fuselage bend- 
ing modes of the unheated vehicle are considered and 
are given by = 2.95 Hz, o >2 = 5.72 Hz, and u >3 = 
7.74 Hz. The damping ratio for structural vibration 
modes is typically small and was chosen to be = 
0.01 for all three frequencies. The frequency response 
of the resulting flexible mode function in Eq. (32) 
is shown in Fig. 5. Due to the “uniform beam” ap- 
proximation an amplification was necessary to achieve 
a reasonable response magnitude which in this case 
translates into 0.25° peak uncertainty in a for a 1° si- 
nusoidal variation in elevon at w = Also given in 
Fig. 5 is the response of the rigid body due to elevon 
deflections and a cover function which was used for 
the controller design. The cover function is given by 


WW = 


*C (* + 2) 

S 2 + 2 Cc^cS + u 2 


(35) 


with k c = 2.5, £c = 0.25 and w c = 19.53 rad/sec. 

The effect of aerodynamic heating during high- 
speed flight through the atmosphere on the elastic 
characteristics of a hypersonic vehicle was addressed 
in Ref. [13]. One conclusion was that for increasing 
aerothermal heating loads (related to increasing Mach 
numbers) the natural frequencies of the fuselage bend- 
ing modes decreased. This envokes the potential of 
rigid/elastic body mode coupling. Although the longi- 
tudinal model used in the present study accounts only 
for one Mach number, this effect will be considered by 
letting the frequencies of the elastic modes vary from 
100% to 60% of their values for the unheated vehicle. 
The cover function used for the controller design was 
fitted accordingly by adjusting the parameters from 
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(dash), and cover function (dash dot) 


Eq. (35). One feature of the elastic mode approxima- 
tion used in this study is the fact that for decreasing 
frequencies the peak value of its frequency response 
increases. Thus, the closer the elastic modes get to 
the rigid body frequencies, the more uncertainty in 
the angle of attack is inserted into the system. This 
uncertainty is labeled “uncertainty 2” in Fig. 4. 

The third type of uncertainty introduced into the 
system is uncertainty on the control effectiveness in 
both elevon and equivalence ratio control channels, 
depicted “uncertainty 3” and “4” in Fig. 4. The un- 
certainties are represented multiplicatively around the 
actuators and occur simultaneously. This issue was 
already addressed in Ref. [8] using constant levels 
of uncertainty, here, the uncertainty weight on each 
channel is modeled as 


__ kgctjs 4 - 10 ) 

s + 1000 


(36) 


to account for increasing uncertainty with increasing 
frequency at the plant input. k aci determines the un- 
certainty level for low frequencies. 


5. Preliminary Sensitivity Studies 


Before analyzing the impact of the various uncer- 
tainties introduced in the previous section on the 
achievable robust performance, an Hoc-controller for 
the nominal model without any uncertainty is de- 
signed in order to scale nominal performance. The 00 - 
norm from disturbances to errors disregarding uncer- 
tainty for the closed loop system was ||T e d||oo = 0.78. 

In a first series of investigations only the effect of 
parameter uncertainty in cm a was addressed. For a 
variety of uncertainty levels an Hoo-^mtroller was de- 
signed and a series of succesive D-K iterations lead to 
the final /i-controller. Fig. 6 shows the curves for the 
upper bound on the p values representing robust per- 
formance bounds of the system for designs with 25%, 



Figure 6: Robust performance bounds for 25% (dot), 
50% (dash dot), 75% (dash), and 100% (solid) uncer- 
tainty in Cm a- 

50%, 75%, and 100% uncertainty in the aerodynamic 
coefficient cm<*- The values of the /i bounds do not 
exceed 1 implying that robust performance is achiev- 
able even for 100% uncertainty in cm«. 

Similar to the cm a uncertainty study, investigations 
with respect to the angle of attack uncertainty due to 
fuselage bending were carried out. As discussed in 
the previous section, the frequency range over which 
the flexible modes occur is varied to account for aero- 
dynamic heating effects. Robust performance bounds 
for angle of attack uncertainties specified in Table 1 
are shown in Fig. 7. Again, robust performance was 
achievable in all cases. However, the results indicate 
that achievable robust performance is very sensitive 
to the frequency of the first flexible mode. 

The last uncertainty representation to be analyzed 
is the simultaneously occuring uncertainty in the actu- 
ators. Low frequency uncertainty levels of 10%, 20% 
and 30% were imposed on the control effectiveness. 
The corresponding robust performance bounds are de- 
picted in Fig. 8. Note that robust performance is not 



Figure 7: Robust performance bounds for uncertainty 
in a: case 1 (dash dot), case 2 (dash), case 3 (solid). 
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case 

frequencies 
(rad /sec) 

% of unheated 
frequencies 

cover function 
*c Ce w c 

i 

18.53, 35.34, 48.63 

100% 

2.5 

0.25 

19.53 

2 

14.82, 28.27, 38.91 

80% 

3 

0.25 

15.82 

3 

11.12,21.20,29.18 

60% 

3 

0.2 

11.12 


Table 1: Uncertainty levels in a. 



Figure 8: Robust performance bounds for 10% (dash 
dot), 20% (dash) and 30% (solid) in the actuators. 

achievable at the 30% level. The simultaneous pres- 
ence of two uncertainties clearly reduces the level of 
obtainable robust performance. However, it should 
be noted that uncertainty in controlling thrust may 
be less than that for controlling pitching moment. 

The sensitivity studies indicate that uncertainty in 
cm<* is not as critical as uncertainty in a or in the 
control effectiveness. The effect of propulsion system 
variations on the pitching moment due to angle of 
attack changes of the rigid body are not potentially 
“dangerous” for the stability of the vehicle. Flexible 
fuselage bending becomes more important the closer 
the elastic frequencies get to the rigid body frequen- 
cies. Uncertainty in the control effectiveness is also 
crucial for higher uncertainty levels at low frequencies. 
These results illustrate the trends when the uncertain 
effects appear one at a time. A design allowing for 
mutual coupling is investigated in the following sec- 
tion. 

6 . Design for Robust Performance 

A ^-controller for all types of uncertainty being 
present simultaneously was designed. The results of 
a design with 25% uncertainty in cm ay uncertainty in 
the angle of attack corresponding to case 1 in Table 1, 
and 10% uncertainty in the actuators are illustrated 
in Fig. 9. Robust performance is achieved at all fre- 
quencies. The time responses for the nominal model 



Figure 9: Robust performance (solid), nominal per- 
formance (dash) and robust stability (dash dot) for a 
fi design including all uncertainties. 

with no uncertainty for simultaneous step commands 
in velocity (V c = 100 ft/sec) and altitude (h c = 1000 
ft) while encountering longitudinal and vertical at- 
mospheric turbulence are shown in Figs. 10-15. Ve- 
locity and altitude responses compare well with the 
Hoo design for the nominal model (Figs. 10 and 11). 
Both meet the performance requirements. The alti- 
tude curve exhibits a little ‘dip’ towards the end which 
is due to influence of turbulence. The magnitudes of 
the variations in angle of attack, normal acceleration, 
elevon deflection and fuel flow are reasonable consid- 
ering the level of turbulence present in the simulation. 
Figs. 12-15 illustrate the effect of turbulence on the 
system. 

Fig. 13 shows that turbulence causes a ±0.05g vari- 
ation in the normal acceleration sensed by the pilot. 
This can be reduced by penalizing n t in the design 
at the expense of increasing the bandwidth of the 
pitch loop and increasing the coupling with the flexi- 
ble mode. 

The singular values of the loop transfer function 
with the loop broken where velocity and altitude are 
fedback to the controller are shown in Fig. 16, for 
both the Hoo and the /i design. The bandwidth of 
the /i-Bynthesis design is reduced and the crossover 
freequencies of the minimum and maximum singular 
values are closer together. 

Figs. 17 and 18 compare the time responses of a 
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flmt (mc) 


Figure 10: Velocity response for nominal perfor- 

mance, Hoo design (solid), y design (dash). 


40 60 

tmc (mc) 


Figure 13: Normal acceleration response for y design, 
with (solid) and without (dash) turbulence. 


ton* (mc) 


Figure 11: Altitude response for nominal perfor- 

mance, Hoo design (solid), y design (dash). 


time (mc) 


Figure 14: Elevon response for y design, with (solid) 
and without (dash) turbulence. 


*ne (mc) 


Figure 12: Angle of attack response for y design, with 
(solid) and without (dash) turbulence. 


40 60 

time (mc) 


Figure 15: Fuel equivalence ratio response for y de- 
sign, with (solid) and without (dash) turbulence. 









Figure 16: Singular value plots for loop broken at V 
and h feedback, Hqq design (dash), /i design (solid). 



Figure 17: Velocity responses for Hqo design (solid) 
and perturbed fi design (dash). 

perturbed system using the ^-controller to the nom- 
inal response of the Hqo design. The perturbations 
that are introduced into the system were chosen to be 
the worst case and of the same magnitude as those 
for which the /j-cont roller was designed. Figs. 17 and 
18 illustrate that performance can be maintained 
by a ^-controller for a system with four simultaneous 
perturbations. 

The effects of the uncertainties cannot be consid- 
ered seperately. Uncertainty in control effectiveness 
changes the elevon deflections affecting the angle of 
attack changes due to flexible body bending. This 
variation in a adds to the rigid body angle of at- 
tack which again influences the impact of the propul- 
sion system on the pitching moment. Obviously, cou- 
pling between the individual effects occurs affecting 
and possibly amplifying the uncertain characteristics. 
When a controller design for higher levels of uncer- 
tainty is performed, robust performance cannot be 
achieved (see Fig. 19). 100% uncertainty in un- 
certainty in a correponding to case 3 in Table 1 and 
20% uncertainty in control effectiveness were used. 



Figure 18: Altitude responses for Hoo design (solid) 
and perturbed u design fdashV 



Figure 19: Robust performance (solid), nominal per- 
formance (dash) and robust stability (dash dot) for a 
/i design with increased uncertainty levels. 

These vadues allowed robust performance in the indi- 
vidual designs. Still, with a maximum value of the 
robust performance bound at around 1.35, the result 
is not “devastating”. Robust stability is maintained. 

7. Controller Order Reduction 

The order of the /i-controller accounting for all un- 
certainties is 38. The vehicle was represented by a five 
state longitudinal model, and the generalized plant 
including all frequency dependent weights for perfor- 
mance and uncertainty is 16 th order. Designing for ro- 
bust performance considering structured uncertainty 
(D scaling) increased the order of the controller con- 
siderably. When implemented, large order controllers 
can create time delays which may be undesirable. One 
solution to this problem is to use model order re- 
duction on the controller realization. This method, 
though, does not consider the properties of the closed 
loop system when reducing the order of the controller, 
i.e. robustness properties are not guaranteed. 

To investigate this issue, an order reduction of the 
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Figure 20: Robust performance bounds for 26 state 
reduced order controllers. 

controller was performed using several different meth- 
ods provided by the MATLAB toolboxes p-TOOLS 
[14] and Robust Control Toolbox [15]. Not all meth- 
ods could be applied since the controller is unstable. 
In a first step, the order of the controller was reduced 
by 12 resulting in a 26 state controller. Robust perfor- 
mance bounds are illustrated in Fig. 20, which clearly 
demonstrates that robust performance of the full or- 
der controller is not maintained. The methods dis- 
played in Fig. 20 are truncating (sirunc) and residual- 
izing (sresii) a controller which has been transformed 
to bidiagonal form (sfrans) Ref. [14]. The other meth- 
ods utilized in this comparison are a balanced model 
reduction ( balmr ) and an optimal Hankel approxima- 
tion (ohkl) from Ref. [15]. Both preserved robust 
performance for the 26 state controller (not shown) 
but lost this property when the controller order was 
chosen to be 20, (see Fig. 21). 

Another approach to design low order compen- 
sators is to constrain the order of the controller in 
the design process, and promising results have been 
achieved in Refs. [16] and [17]. This technique com- 
bined with the framework of p synthesis will be ad- 
dressed in future research. 

8. Conclusions 

A control study of an hypersonic vehicle with air- 
breathing propulsion has been performed in this pa- 
per. Special consideration is given to modeling the un- 
certainty due to propulsive and aeroelastic effects and 
control effectiveness. The impact on achievable ro- 
bust performance was examined. A sensitivity study 
of robust performance achievable with various levels of 
uncertainty shows that flexible fuselage bending and 
uncertainty in control effectiveness have a greater im- 
pact on achievable robust performance than propul- 
sion system variations due to angle of attack changes. 
A design example for robust performance illustrated 



Figure 21: Robust performance bounds for 20 state 
reduced order controllers. 

that Hoo and ^i-synthesis is ideal for addressing the 
problem of multiple uncertainty sources present in 
flight control systems. Acceptable uncertainty levels 
degrade when all uncertainties occur simultaneously. 
This implies that uncertainty modeling will play a cru- 
cial role in successful flight control system design for 
hypersonic vehicles. 

An order reduction of the resulting controller dis- 
closes that robust performance is not guaranteed for 
the reduced order system. Future research will focus 
on fixed order controller design and the issue of treat- 
ing real parameter uncertainty. 
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The effect of turbulence on the system is modeled 


by 


where 


x = Ax + Bu + B g 
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The elements in B g are the corresponding elements of 
the A matrix. The pitch rate due to gust q g can be 
expressed in terms of the angle of attack rate 


w. 


The numerical values for the A and B matrices are 
given by 
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